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Abstract. Given a multivariate generating function F(zi, . . . , Zd) = X! a ri,— .^Zi 1 ■ ■ ■ z^ d , we determine 
asymptotics for the coefficients. Our approach is to use Cauchy's integral formula near singular points of 
F, resulting in a tractable oscillating integral. This paper treats the case where the singular point of F 
is a smooth point of a surface of poles. Companion papers will treat singular points of F where the local 
geometry is more complicated, and for which other methods of analysis are not known. 



1. Introduction 

The generating function F(z) := Ylr^=o arzT f° r the sequence 00,0,1,02,... is one of the most useful 
constructions in combinatorics. If the function F has a simple description, it is usually not too hard to 
obtain F as a formal power series once one understands a recursive or combinatorial description of the 
numbers {a r }. One may then analyze the analytic properties of F in order to obtain asymptotic informa- 
tion about the sequence {a r }. While still part art and part science, this latter analytic step has become 
quite systematized. Stanley! (|1997| ) in his introduction to enumerative combinatorics gives the example 

2 

of the function F(z) = exp(z + from which he says "it is routine (for someone sufficiently versed in 
complex variable theory) to obtain the asymptotic formula a r = 2~ 1 / 2 W 2 e~ r / 2+v/ ^~ 1 / 4 ." Routine, in this 
case, means a single application of the saddle point method. When F has singularities in the complex 
plane, the analysis is often more direct: the location of the singularities and the behavior of F near these 
determine almost algorithmically the asymptotic behavior of the sequence {a r }. For those not sufficiently 
versed in complex variable theory, two useful sources are Henrici| ( |19771 ) and pdlyzko (1995). The transfer 
theorems of Flajolet fc Odlyzkc] ( |1990[ ) encapsulate much of this knowledge in a very useful way; see also 
Will (1994) for an elementary introduction. 

When the sequence 00,01,02,... is replaced by a multidimensional array {a rij ... jrd }, things become 
much more hit and miss. Let us use boldface to denote vectors in C d or N d , and use multi-index notation, 



■ I'd 



and z r denotes the product z x 



so that a r denotes the multi-index a 
write in expanded form for clarity. The generating function F : 
one-dimensional generating function by 



• • • z. 



which we will sometimes 
C is defined analogously to the 



F(z) 



2 ar 

r£N d 



Surprisingly, techniques for extracting asymptotics of {a r } from the analytic properties of F were, until 
recently, almost entirely missing. In a survey of asymptotic methods, Bender (1974) says: 
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Practically nothing is known about asymptotics for recursions in two variables even when a 
generating function is available. Techniques for obtaining asymptotics from bivariate generating 
functions would be quite useful. 

In the intervening 25 years, some results have appeared, addressing chiefly the case where the array {a r } 
obeys a central limit theorem. Common to all of these is the following method. Treat {a r } as a sequence 
of (d — l)-dimensional arrays indexed by r^; show that the n th (d — l)-dimensional generating function 
is roughly the n th power of a given function; use this approximation to invert the characteristic function 
and obtain a Central Limit Theorem. We refer to these methods as GF-sequence methods. The other 
body of work on multivariate sequences, which we will call the diagonal method, is based on algebraic 
extraction of the diagonal, as found in Hautus fc Klarner ( 1971| ) ( see also Furstenberg ( 1967 ) and later 
[Lipshitzj ( 1988|) for an algebraic description of the scope of this method; variants are described in Stanley 
(|1999| ) and TPippengerj (^500|) ). 

The most fundamental GF-sequence result is probably Bender &: Richmond! ( |1983j ), with extensions 
appearing in later work of the same authors. [Flajolet fc Sedgewick ( 1997| ) present a version of the 
same idea which holds in much greater generality. Gao & Richmond ( 1992| ) go beyond the central limit 
case, using the transfer theorems of |Flajolet fc Odlyzko ( 1990 ) to handle functions that are products of 
powers with powers of logs. Recent work of Bender and Richmond (Bender & Richmond 1996, Bender &| 
Ric hmond 1999| ) extends the applicability of the central limit results to many problems of combinatorial 
interest; see also ( |Hwang 1995| , |Hwang 1998 where more precise asymptotics are given, and Hwang] 
( |1998o| ), which extends some results to the combinatorial schemes of [Flajolet &: Soria] ( |1993j ) . This does 
not exhaust the recent work on the problem of multivariable coefficient extraction, but does circumscribe 
it. 

The present paper, together with forthcoming companion papers, takes aim at a large class of multivari- 
able coefficient extraction problems, for which a fair amount of information can be read off in a systematic 
way. An ultimate goal (not our only goal) is to systematize the extraction of multivariate asymptotics 
sufficiently that it may be automated, say in Maple. Everything we do, we do with complex contour 
integration. In this regard, our methods are most similar to those of Bertozzi &: McKenna ( 1993] ), who, 
as we do, provide a general framework for harnessing the multivariable theory of residues for exact and 
series computation of coefficients. A more detailed description of our method will be given in Section ||, 
but here is an outline. 



(1) Use the multidimensional Cauchy integral formula to represent a r as an integral over a 
d-dimensional torus inside C d . 

(2) Expand the surface of integration across a point z where F is singular, and use the 
residue theorem to represent a r as a (d — l)-dimensional integral of one-variable residues. The 
choice of z determines the directions in which asymptotics may be computed. 

(3) Put this in the form of an integral J exp(A/(z))^(z) dz for which the large-A asymptotics 
can be read off from the theory of oscillating integrals. 

In the rest of this introductory section, we describe the scope of our methods. Figure |l| depicts a 
classification of generating functions and illustrates the remainder of this paragraph. If a formal power 
series is nowhere convergent, analytic methods are useless. Among those power series converging in some 
neighborhood of the origin, there are three possibilities: a function may be entire, may have singularities 
around which analytic continuations exist, or it may be defined only on some bounded subset of C d . Our 
methods are tailored to the second class. The third class, although in some sense generic, seldom arises 
in any problem for which the generating function may be effectively described. Incomplete asymptotic 
information is available via Darboux' method; details of this method in the univariate case are given in 
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Figure 1 . Classification of generating functions 



Henrici| (|1977|) and |Odlyzkoj QL99SQ . The first class can and does arise frequently. Our methods are simply 



not equipped to handle entire functions, and systematizing the asymptotic analysis of coefficients of entire 
generating functions remains an important open problem. 

For the remainder of this paper, we will assume that the formal power series F converges in a neigh- 
borhood of the origin and may be analytically continued everywhere except a set V of complex dimension 
d — 1 which we call the singular variety. The point z in step 2 is an element of V, and the behavior of 
V near z greatly affects the subsequent analysis in step 3. This paper addresses the case where z is a 
smooth point of V at which F has a pole. The forthcoming companion papers will address cases where z 
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is a multiple point or a cone point. We do not know whether cases where z is a cusp of V arise, but if so, 
the subsequent analysis has mostly been carried out in the work of AnioFd, Guscm-Zadc fc Varchenkcj 

( HP - 

The chief purpose of this study is to give a solution to the problem of asymptotic evaluation of coef- 
ficients that is as general as possible. An important part of this is re-derivation in a general setting of 
results obtainable via GF-sequence or ad hoc methods. We show in Section || how unifying these results 
allows us to show that our method successfully finds asymptotics for every function in a certain large class. 
Familiar examples from this class include: lattice path counting, various known generating functions for 
polyominos and stacked balls, enumeration of Catalan trees by number of components or surjections by 
image cardinality (see Flajolet & Sedgewick ( 1997| ), stopping times for certain random walks (see Larsen 
& Lyons (|199£ )), as well as the examples given in the GF-sequence papers of Bender ( 1973| ) and |Bendci 
&; lichmond| ( 1983j ): ordered set partitions enumerated by number of blocks, permutations enumeratec 
by rises, and Tutte polynomials of recursive sequences of graphs. 

Nevertheless, our pursuit of this problem was also motivated by some specific applications which we 
mention briefly now and discuss more thoroughly later. These are cases where known methods do not 
suffice to obtain complete asymptotic information. There is a class of tiling enumeration problems for 
which an explicit three variable rational generating function may be obtained. This class includes the 
Aztec Diamond domino tilings of Cohn, Elkies fc Propp| ( |1996| ). Asymptotics in the so-called region of 
fixation are obtained from analysis of the smooth points of V (Theorem |3.5| below), while asymptotics 
in the region of positive entropy are derived from analysis of the cone point. |Cohnfc Pemantle| (|2000|) 
applies a cone point analysis to a tiling enumeration problem for which the only previous results are some 
pictures via simulation ( http: / / www.math.harvard.edu / ~cohn /picture.gil ) . Another motivation has been 
to solve the general multivariable linear recursion. Depending on whether one allows forward recursion in 
some of the variables, one obtains either rational or algebraic generating functions. The general rational 
function may have any of the types of singularities mentioned above: smooth points, nodes, cones, cusps, 
branchpoints, etc. Even the simple rational generating function 1/(3 — 3z — w + z 2 ) of Example |3.4| req uires 
two separate analyses in order to get asymptotics in all directions. We will see that Theorem |3.1| gives 
asymptotics in one region, while Theorem 3.3 is required for other directions. 

Asymptotics derived near smooth pole points nearly always exhibit central limit behavior. Smooth 
pole points are the topic of this first paper, and are exactly the case to which existing methods may 
apply. While one function of this paper is to lay foundations for the cases in which the singularity is more 
complicated, there are several ways in which it improves upon available analyses of the smooth case. 

First, most of the existing results assume that the singular point z G V has positive real coordinates, 
and that it is strictly minimal in a sense defined in the next section. This assumption often holds when 
the coefficients {a r } & r e nonnegative reals, though it will fail if, for example, there is any periodicity. The 
assumption always fails when the coefficients {a r } have mixed signs, as is the case for example with the 
generating functions (1 — zw)/(l — 2zw + w 2 ) and 1/(1 — 2zw + w 2 ) for the Chebyshev polynomials of 
the first and second kinds (Comtet 1974 , page 50). GF-sequence methods may be adapted to some of 
these situations. Indeed, the presentation of these methods by Flajolet &: Sedgewick ( 1997 , Theorem 9.7) 
accomplishes this adaptation in great generality. But certainly there are cases such as the rational 
generating function 1/(1 — z — w + j3zw), where the points z with given moduli form a continuum and 
standard GF-sequence methods are not sufficient. 

Second, our methods obtain automatically a full asymptotic expansion of a rii ___ jTd in decreasing powers 
of the indices rj. This is certainly not inherent in the existing results, whose relatively short proofs involve 
inversion of the characteristic function (see however Hwang! ( |1995| ) and [Hwang! ( 1996| ) for something in 
this direction). The expansion to n terms is completely effective in terms of the first n partial derivatives 
of 1/F at z, as is the error bound. 
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Third, these results explicitly cover the case where the pole at z has order greater than 1. The behavior 
in this case is not according to the central limit theorem. The only existing work addressing this case is 
pao fc Richmond] ( |1992 ), and they require nonnegativity assumptions, as mentioned above. In the case 
where F = G k is an exact power, one could attempt first to solve the problem for G and then to take the 
A;-fold convolution. This is much harder than the present approach, as may be seen by the rather involved 
computation in Cohn et al. (1996). 

Fourth, the potential for increasing the scope to new applications seems greater for contour methods 
than for GF-sequence methods. The contour method reduces the asymptotic problem to the problem of an 
oscillating integral near a singularity, which can almost certainly be done. By contrast, the GF-sequence 
method requires first an understanding of the sequence of (d — l)-dimensional generating functions aris- 
ing from the given (i-dimensional generating function, and then another result in order to transfer this 
information to asymptotics of the coefficients a r . 

Fifth, although our results in the case of smooth pole points are often similar to those obtained by 
GF-sequence methods, our hypotheses are quite different. In Section || we show how our hypotheses may 
be universally established for functions that generate nonnegative values and are meromorphic through 
their domain of convergence. 

Finally, we compare our method to recent results from the diagonal method. It is known ( Lipshitz 1988 ) 
that the diagonal sequence a n ,n,...,n of a multivariate sequence with rational generating function has a 
generating function satisfying a linear differential equation over rational functions. Much is known about 
how to compute this equation (see for example phyzak Salyy| ( |1998[) ). If one wants asymptotics on the 
diagonal, or in any direction where the coordinate ratios are rational numbers with small denominators, 
then these methods give results that are in theory at least as good as ours. The method, however, is 
inherently non-uniform in the direction, so there is no hope of extending it to larger sets of directions, 
which is what we accomplish in the present work. 

The remainder of the paper is organized as follows. In the next section we set forth notation and define 
the terms necessary to state the main results of the paper. The main results are stated in Section ||, 
and examples are given. The next section contains a proof of these results, modulo the computation of 
some oscillating integrals. This computation is carried out in Section [|. Section [6| outlines some details 
of taxonomy and discusses universality of the method of complex contour integration. The final section 
states some open problems. 



2. Notation and Preliminaries 

The main results of this paper give asymptotics valid under certain geometric assumptions on V and 
computable from some quantities that are in turn effectively computable from the generating function F. 
Thus in addition to setting out basic notation, we need to define some terms related to the geometry of 
V and some quantities associated with F. 

2.1. Notation. Throughout the paper, F = ^a r z r will denote a function on C d analytic in a neighbor- 
hood of the origin. The (open) domain of convergence of the power series will be denoted T>. For z £ C d , 
let T(z) denote the torus consisting of points w with \wj\ = \zj\ for 1 < j < d and let D(z) denote the 
closed poly disk of points w with \ wj\ < \zj \ for 1 < j < d. Recall (see Hormander] (|1990| )) that the domain 



T> is a union of tori T(z) and is logarithmically convex, that is, the set 

log£> := {xel d : (e Xl , e Xd ) £ V} 

is a convex subset of R d and is an order ideal, that is, it is closed under < in the coordinatewise partial 
order. 
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We assume throughout that F = G/H, where both G and H are analytic in a neighborhood of D(z) 
for some point z. In particular, every meromorphic function satisfies this condition The set where H 
vanishes will be denoted V. Many of our examples will be in dimension 2, in which case we will often use 
z and w in place of z\ and Z2, use (z,w) in place of z, and use (r,s) in place of (n,^). We sometimes 
need to treat C d as C d_1 x C (although symmetry of the coordinates is preserved most of the time). 
Accordingly, when the dimension is greater than 2, we use z to denote (z%, . . . , z^—i)- Partial derivatives 
will be denoted Hi for |^ and so forth; in dimension 2 we will also use H z and H w . 

As is usual for asymptotic analyses, we let / ~ g denote f/g — > 1, with the limit taken at infinity 
unless otherwise specified. The function / is said to be rapidly decreasing if f(x) = 0(x~ N ) for every N, 
and is said to be exponentially decreasing if f(x) = 0(e~ cx ) for some c > 0. We also use the symbol "~" 
to denote asymptotic expansion. Thus 

/ ~ hn 9n 

is normally taken to mean that / — 5^n=0 &nSn = (°n9n), where b n € C and {g n } is a fixed sequence of 
functions such that g n +\ = o{g n ) for each n. We broaden this to allow b n = when n ^ 0, so that the 
remainder term need only be o(g n ) and not o(b n g n ). In particular, if 

oo 

f( X ) ~g( X ) .J2c n X- n 
n=0 

with Co = 1, then we say we have obtained a full asymptotic expansion for / in decreasing powers of x 
with leading term g. 

2.2. Geometry of V. As in the one-dimensional case, the points of V nearest the origin are the most 
important. Accordingly we define a point z G V to be minimal if Vf)D(z) C T(z); we say that z is locally 
minimal if the analogous relation holds with V replaced by a neighborhood of z in V. Divide the minimal 
points of V into three types. Say that z is strictly minimal, finitely minimal or toral, according to whether 
the cardinality of V Pi D(z) is 1, finite, or infinite. When infinite, the intersection must be uncountable. 
If z is a minimal point of V then the interior of D(z) is contained in D, so the assumption that G and H 
are analytic on a neighborhood of D(z) is just a little stronger than what is true automatically. 

A simple pole of F is a point z € V where H vanishes to order 1. Equivalently, the gradient V H does 
not vanish. Let z be a simple pole of F and assume for specificity that is nonzero at z. By the implicit 
function theorem, there is a neighborhood of z where V may be parametrized by z^ = g{z\, . . . , z^-i) for 
some analytic function g. We will always use g to denote this parametrization. 

We will see later (in the proof of Theorem |6.3| ) that under some hypotheses on F, minimal points of V 
are always found in the positive real orthant. A relation true in complete generality is the following. 

Lemma 2.1. Let z be a simple pole of F and suppose that ZdH^ does not vanish there. If z is locally 
minimal then for all j < d, the quantity ZjHj/(zdHd) is real and nonnegative. 

Proof. Given 9 and j, let z^ be the result of varying z by multiplying the j th coordinate by e lS and 
adjusting the last coordinate so as to remain on V (that is, 

zP = g(z\,... , Zj-i, Zje lB , zj + i, . . . ,Zd-i))- Differentiating the relation H(z^) = implicitly with 



1 The greater generality allows us to cover examples such as th e generating function for self-avoidin g random walks ( Chaves 



ic Chayes 1986| ) or percolation paths in the subcritical regime (Campanino, Chayes & Chaycs 1991). In these cases, all the 



work is in showing the function is meromorphic in a neighborhood of D(z). Without further knowledge, the authors then 
conclude central limit behavior. 
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respect to 9 at yields 

dz {9) 

(2.1) izjHj + Hd ^ = _ 

(0) (0) 

By minimality of z, we know that the modulus of z d has a minimum at 9 = 0, hence (dz d /d9)/zd is 
purely imaginary. Plugging this into (2J.) proves that ZjHj /(zdHd) is real. If ZjHj/(z c iH c i) = —j3 < 



then V has a tangent vector at z in the direction —zjej — flzjed, where ej is the j th coordinate vector. 
This contradicts minimality. Hence ZjHj/ (zdHd) > 0. □ 

Definition 2.2. Define dir(z) to be the equivalence class of (complex) scalar multiples of the vector 
(z\Hi, . . . , ZdHd), defined whenever ZjHj does not vanish for all j. By the previous lemma, when z is a 
minimal pole of F with nonzero coordinates, dir(z) is a well defined element of MP d_1 . 

The importance of dir is that analysis of F near z yields asymptotic information about a r with r £ 
dir(z). The function dir appears in GF-sequence method literature as m. When z G &D is on the 
boundary of the domain of convergence, dir(z) is the normal to the support hyperplane of the convex set 
\ogV at the point (log \z\\, . . . ,log|^|)- 

We now define a few more quantities associated with F and g. Again, we will reserve the names of 
these functions, so as not to burden the notation with subscripts and arguments. If z is a simple pole of 
F with ZdHd not vanishing there, define a function ip on a neighborhood of z by 

(2.2) ^(w) = - lim (w - g(w)) F(yr,W \ 

Suppose now that w € T(z) and write Wj = Zje 1 ® 3 . For fixed r with ^ 0, define a function / on a 
neighborhood of z in T(z) by 



(2.3) /( * )=log (^|l) +iE 



-1 



We will be parametrizing integrals over T(z) by 9, so we will want the above function expressed 
in terms of 9. We therefore compose with the map M taking 9 to w defined by M(9\,... ,9d-i) = 
(zie tdl , . . . , Zd-ie 1 ^' 1 ), and define the functions g := g o M, / := / o M, tp := ip o M . 



Although it is not obvious yet, / will always vanish at to at least two orders (Lemma 4.2 below), and 



the hypothesis Q in Theorem 3J is equivalent to / having nonvanishing quadratic term. For ease of 
reference, Table |2.2| summarizes the foregoing definitions, stratified by how many times the given data G 
and H have been manipulated. 

3. Statement of results, with examples 

Before going on, we pause to state a prototype of our results in the simplest possible setting, namely 
where the number of variables is 2, the functions G and g are as nondegenerate as possible, and only the 
leading term asymptotic is given. The proof is in Section ||. 

Theorem 3.1. Let F = G/H be a meromorphic function of two variables, not singular at the origin. 
Define 

Q(z, w) := -w 2 H 2 w zH z - wH w z 2 H 2 z - w 2 z 2 {H 2 H ZZ + H 2 H WW - 2H Z H W H ZW ) . 

Then 

G(z,w) _ r _ s -wH w 
■ -z w 



V2tt ' V S Q 



8 



ROBIN PEMANTLE AND MARK C. WILSON 



Table 1. Reserved notation in remainder of this article 



Given information: 

the function F in the form G/H 



First level: 

g parametrizes the zero set, V of H 

dir(z) is the coordinatewise product (V H) ■ (z) in projective space 

Second level: 

ip is the residue in Zd of F/zd at points (z,g(z)) 

f is log g, plus a term linear in log Zj and depending on r. 

Third level: 

ip,g and / are ip,g and / expressed in terms of 9 



uniformly as (z, w) varies over a compact set of strictly minimal, simple poles of F on which Q and G 
are nonvanishing, and (r,s) G dir(z,w). 

Remarks: Usually the expression in the radical will be positive real, as will the coefficients a rs . The result 
is true in general, though, as long as the square root is taken to be —wH w times the principal root of 
Q/(—wH^). Also note that when (r, s) G dir(z,w) then the expression wH w /s is coordinate-invariant, 
that is, equal to zH z /r. Thus the given expression for a r>s has the expected symmetry. 

Example 3.2 (Lattice paths). 

Let a T)S be the number of nearest-neighbor paths from the origin to (r, s) moving only north, east 
and northeast; these are sometimes called Delannoy numbers Stanley (|1999| , page 185). The generating 



function is F{z, w) = l/(l—z—w—zw). The zero set V of H = 1 — z—w—zw is given by w = (l — z)/(l+z) 
and the minimal points of V are those where w G [0, 1]. With the help of relations that hold when z G V 
we may compute as follows. 

H z = -1-w 
—zH z = 1 — w 

Q = (1 - z)(l - w)(l - zw) 
zH z 1 — w 1 — w 2 



wH w 1 — z 2w 

with H w and —wH w given by reversing z and w. Ks z varies over [e, 1 — e], the functions Q and G := 1 

do not vanish. The minimal pair (z,w) that solves (r,s) G dir(z,w) is given by z = (vV 2 + s 2 — s)/r and 

w = (\/r 2 + s 2 — r)/s. Theorem |3.1| then gives 

/ \ — r / \ —s i — 

I \fr 2 + s 2 — s\ I \r 2 + s 2 — r \ 1 



2ir V s 1 — zw 



y/r 2 + s 2 — s \ I \Jr 2 + s 2 — r\ /l 



rs 



2vr Y ( r + s - Vr 2 + s 2 ) 2 Vr 2 + s 2 
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uniformly when r/s and s/r remain bounded. In particular, when r = s = n, this gives the following 
formula for the n th diagonal coefficient (which may alternatively be obtained by computing the diagonal 
generating function (1 — 6s + s 2 )~ 1 / 2 according to the method given in |Stanley| ( jjjgg , Section 6.3): 

T 2-v 4 

2tt2 - y/2 




The computations in Theorem in terms of the values and derivatives of G and H are explicit. As 
we state more general theorems, it becomes cumbersome and in fact obfuscating to give formulae for the 
expansion coefficients directly in terms of derivatives of G and H. This is one reason we have already 



introduced the functions in Table 2.2. It should be emphasized, however, that while we use higher level 
quantities in the statements of subsequent theorems, each expansion coefficient can be computed from 
finitely many derivatives of G and H. We begin with a relatively explicit computation for the general 
two- variable case. 

For k at least 2, we define constants 

(3.1) A+(k,l) := \r n + 1 

(3.2) A(k,l) := ir 0-^) ( 1 + e^A^W-^) ) if fc is odd, 

(3.3) A(k,l) := — r ( — j- — ) if k, I are even, 



Let 



k 

A(k,l) := if A; is even and I is odd. 

*.) = /(»)"» ^.(i + M^)*, 

where c& is the first nonvanishing Taylor coefficient of f{x) = YlJ^k c j x ^ anc ^ the argument of c]/ k is 
taken between —ir/(2k) and n/(2k). Let rj denote the inverse function to y and let {bj} be the Taylor 
coefficients of (ip o rj) ■ rj 1 . Clearly each {bj} is determined by finitely many partial derivatives of G and 
H , and the index Iq of the first nonvanishing b\ is the same as the order of vanishing of tp at 0. The 
coefficients bi are easily computed from the coefficients bj := i/!«'(0)/j! and Cj := f^'{0)/jl; in particular, 
if / ~ CkX k near then 



-l/fc 



(3-4) 6 Jo = b h) c k . 

Theorem 3.3. Let F = G/H = ^a rs z r w s have a strictly minimal, simple pole at (z,w). Let k be the 
order of vanishing of f at 0. Let Iq be the order to which G vanishes near (z, w) on V, that is, the largest 
I such that G(z',w') = 0(\z — z'y + \w — w'\ ) as (z',w') — > (z,w) in V. Then there is a full asymptotic 
expansion 

(3.5) a r , s ~ ^-z~ r w- s V A(k, l)b lS ^ l+1 ^ k , 

l>lo 



where A(k, I) denotes A{k, I) if Imjcfc} > and A(k, I) otherwise. The expansion is uniform as (z, w) 
varies over a compact set of strictly minimal poles with (r, s) G dir(z, w) and k and Iq not changing. 
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Figure 2. V for Example 3.4 



Example 3.4 (Cube root asymptotics). Let F(z,w) = 1/(3 — 3z — w + z 2 ). The set V is the set {w 



3z + 3} and g(z) 



3z + 3. The point (1, 1) is in V, indicating that the maximal exponential 



growth rate will be zero. Indeed, for directions above the diagonal, Theorem 3.1 or |3.3| may be used at 
the minimal points {(z,g(z)) : < z < 1}, while each direction below the diagonal corresponds to a pair 
of complex minimal points fitting the hypotheses of Corollary 3/7; the result is that the coefficients decay 
exponentially at a rate that is uniform over compact subsets of directions not containing the diagonal. 

The interesting behavior is near the diagonal. The relevant minimal point is (1, 1), where z r w s = 1 
and the decay is sub-exponential. Computing f"(0) via equation fl4.8|) below gives 



t z(z 2 -4z + 3) 
V -3z + 3) 2 



This vanishes when z = 1, and computing further, we find that / vanishes to order exactly 3 here, with 
C3 := /'"(0)/3! = i. Along with ip(0) = 1, this then results in an asymptotic expansion whose leading 
term is given by 

a rr „ ^(3,0)^/3(1 + e -W3 )r -V3 = E(^3) r -i/3 

2-7T 6V37T 

In Section |7| we discuss the question of computing asymptotics "in the gaps" so as to be able to conclude 



that limsup log a r / log |r| = —1/3 or even limsup|r 



1/3, 



r(2/3) 

6V3tt 



For more than two variables a result holds similar to the two- variable result. 
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Theorem 3.5. Let F = G/H = ^a r z r have a strictly minimal, simple pole at z. Suppose z^Hd does 
not vanish. If the Hessian of f at z is nonsingular, then there is an expansion 



a r ~ z ^ ^ 

where Iq is the degree to which G vanishes on V near the point z . When G does not vanish at z then 
Iq = and 

C = (2^)( 1 - d )/ 2 W- 1 /2^1 



where TL is the determinant of the Hessian at z. 



Example 3.6 (Domino tilings). Random perfect tilings of planar regions by dominos have been a subject 
of some interest, since the analysis by Fisher ( |1961 ) of this model for dimer packing uncovered an exact 
expression for the partition function of the ensemble. A generating function is given in Cohn et al. (1996) 
which allowed the authors to determine, after some cumbersome analysis, which parts of a diamond- 
shaped region (a union of lattice squares approximating the region |x| + \y\ < k) were asymptotically 
deterministic and which contained randomness in the limit as the edge size of the diamond grew. 



An easier analysis in the region of non-randomness is available via Theorem |3.5| together with a slightly 
more informative generating function than was used by Cohn et al.| ( |1996| ). In particular, let 



F(x,y,z) 



oo 

E 

t=0 |r|+|s|<t 



E 



be the generating function for the probability a r>Si t that the tile covering position (r, s) of a random 
diamond of size t will be horizontal. For brevity, we omit formal descriptions of the diamond and its 
indexing. We remark that the use of negative indices (for each fixed t, the sum X)|r|+|s|<< ^ s a polynomial 
in x,x _1 ,y and y~ l ) does not require any alterations in the theory (see Cohn fc Pemantk| ( 2000| ) for 
justification), and that the natural way to parametrize directions is by the pair (r/t,s/t) which varies 
over the diamond \r/t\ + \s/t\ = 1. From Cohn et al. (1996) or from the generation algorithm in Gessel 
[onescu fc Propp ( |1995[ ), one finds 



F(x,y,z) 



z/2 



1 - (x + x~ l + y + y~ 1 )z/2 + z 2 
|Cohn fc Pemantk (2000) show that whenever (r,s,t) satisfy 



t 



s 2 + 2\fr' 1 



+ Ws 2 + 1 



then there is a smooth minimal point (x, y, z) on the pole manifold of F for which (r, s, t) S dir(j;, y, z), 
yielding exponential decay in the direction (r, s, t). The set of directions so parametrized turns out to be 
the region between the diamond \r/t\ + \s/t\ = 1 and the inscribed circle (r/t) 2 + (s/t) 2 = 1/2. Thus they 
recover the description of the region of non-randomness as the complement of the inscribed circle. They 
also obtain descriptions of the region of fixation for related tiling problems in which no other analysis has 
been carried out. 



The extension of all of the above results to finitely minimal points is routine. 
Corollary 3.7. Suppose z is a finitely minimal point ofV with V T(z) = {z±, 

n 
3=1 



,z n }. Then 
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where Ej(r) is the asymptotic expression given by the previous theorems with z = zj. In other words, if 
there are finitely many points on VflT(z), then sum the contributions as if each were strictly minimal. □ 

Example 3.8 (Chebyshev polynomials). Let F(z, w) = 1/(1 — 2zw + w 2 ) be the generating function for 
Chebyshev polynomials of the second kind (Comtet 1974]) ; of course asymptotics for these are well known 



and easy to derive by other means. To use Corollary |3.7| , first find the minimal points for the direction 
(r, s), which are {i{(3 — (3~ l ) /2,ij3) for (3 = =t -y/fipf • Computing Q = 4a 2 (l — a 2 ) and summing the two 
contributions then gives 



r(s — r) 

when r + s is even and zero otherwise, uniformly as r/s varies over compact subsets of (0, 1). 

4. Proofs of main results 
Half of each theorem is easy and follows directly from Cauchy's formula 

1 \ d r 

W J 

it 

where the multi-exponent r — 1 means (7*1 — 1, . . . , — 1). Indeed, if z is a minimal point of V then letting 
T approach T(z) from the inside, we see that |z r |a r does not increase exponentially. If, furthermore, the 
hyperplane through (log \z%\, . . . ,log \zd\) normal to r is not a support hyperplane for logP, then some 
x G \ogD has x • r > (log \z\\, . . . ,log\zd\) • r, and integrating on the torus T(e x ) shows that \z r \a r 
decreases exponentially. All the work, therefore, is in showing the converse, namely that when the 
hyperplane normal to r is a support hyperplane, then z _r does give the right exponential order for o r . 
This is done by evaluating a r . 



(4.1) a r = ( — : ) / w- r - i F(w)dw 



Theorems |3.l| - |3.5| all begin with the reduction of an iterated Cauchy integral to an oscillating integral 



in one fewer dimension. 

Lemma 4.1. Let z be a strictly minimal simple pole of F = G/H . Assume that z^H^ 7^ 0. For a 
neighborhood M of '0 in W 1 ^ 1 define a quantity 

(4.2) S := (2^) 1 - d z- r / exp(-r d /(W(0) dO. 

Then the quantity 



decreases exponentially as M remains fixed and 



00. 



Proof. For e G (0, |^|), let T be the torus T(z) shrunk in the last coordinate by e, that is, the set of w 
for which \uij\ = \zj\, j < d and \ wd\ = \zd\ — s. Write Cauchy's formula as an iterated integral 

(4.3) a r = (-M / w 



V 2vn 



-r-l 



T(z) 



-r d j?, x dw d 



dw . 



Here C\ is the circle of radius \zd\ — e. Let K C T(z) be a compact set not containing z. For each 
fixed w G K, the function F(w, ■) has radius of convergence greater than \zd\- Hence the inner integral 
in equation ( [4.3; ) is 0(|^| + 5)~ Td for some 5 > 0. By continuity of the radius of convergence, we may 
integrate over K to see that 

|z r | [ w- r - 1 F(w)dw 

JKxd 
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decreases exponentially. Thus if M is any neighborhood of z in T(z), the quantity 



-V 

2m J 



w 



F(w) 
¥7+1 



decreases exponentially. Thus we have reduced the problem to an integral over a neighborhood of z. 

Near z there is a parametrization Wd = g(w) of V. Let Ci be the circle of radius \z&\ +e. Then when M 
is sufficiently small compared to e, the image of M under g is disjoint from 62- Fix such a neighborhood. 
For any w G A/", the function F(w, •) has a single simple pole in the annulus bounded by C\ and C2, 
occurring at <?(w). The residue in the last variable of F at g(w) is equal to 



(4.4) 

where ip is defined in (| 



R(w) := —ip(w)g(jw)~ 
2|). Therefore, for each fixed w G J\f, 



F(w 



'V- 



F(w) 

c 2 w d * 



dwd — 2iriR(w). 



But |z r J c _ F(w)dtt; ( i/w r+1 | is bounded by a constant multiple of (1 + e/|2; d |) r ' d (the constant depending 



2 

on the maximum of F on C2) and hence |z r 



(4.5) 



X\ is exponentially decreasing, where 



X = {2my- d / (w)- r_1 5(w)- r ^(w)dw 

= 2vri 1 d z r / — r - ^4 V(w 



Changing variables to Wj = Zje ldj and dvjj = iwjd8j turns the quantity X into 



d-l 



(2vr) 1 - fl! z- r / Y[ e - ir ^i>{e) 



f(w) 



-'V 



d0 



and plugging in the definitions of / and / at fl2,3|) above yields 

(2vr) 1 - d z- r f exp{-r d f(9))tP(0) <W 

which is none other than H. 



AT 



□ 



Remark. It is possible to compute from Cauchy's integral formula in a more coordinate-free way as follows. 
There is a unique holomorphic (d — l)-form wj? on V for which u A dH = G dz\ A • • • A dzj,- Let O be a 
(d+ l)-manifold that is a homotopy from a small torus to a torus at infinity. Then M := On V is a (d— 1)- 
manifold and a r = (27ri)~ d f M w ~ r ~ 1 d-F in the sense of currents, which is none other than f M w r_1 tJi?. 
See Kenyon fc Pemantle| ( |2000[ ) for a more thorough discussion of the foregoing. The manifold M is 



any member of a certain homology class in V with the coordinate axes removed, and choosing M to 
pass through the stationary phase point for the integrand replicates the selection of z with r G dir(z). 
Although more canonical, the coordinate-free method is less suitable for explicit computation, so we do 



not pursue it further here. Suffice it to point out that the conclusion of Theorem 3.5 may of course be 



written in terms more evidently symmetric, as was done in Theorem 3.1 



Equation ( |4.2| ) is easily recognized as the standard form for an oscillating integral. The only unusual 
feature is that the phase is neither real nor purely imaginary. This presents no difficulties, but it does 
necessitate the statement of a result in Section |5| that is a little different from the usual results on purely 
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oscillating integrals, found in, for example, |5tein| ( |1993 ) or Bleistein & Handelsman ( 1986j ). We first 
establish that 9 = is a stationary phase point for the function / when r £ dir(z). 



Lemma 4.2. The quantity /(0) always vanishes. If r € dir(z) then V/(0) 
/ias a strict minimum at 0. 



and irae rea/ port o/ / 



Proof. The first statement is immediate. To prove the second, let j < d — 1 and see from the definition 
of / that 



r dfj(z) 



r d 9j( 



g[z) Zj 

By definition of dir, the ratio rj/(zjHj) is some constant c independent of j, hence 

c-Vd/Cz) = g J {z)H d {z)+H j {z). 

The right hand side of this is the derivative of H(w\ , . . . , Wd—i, s(w)) with respect to at z. By definition 
of g this vanishes, and hence fj(z) = 0. But fj(0) = izjfj(z), so the gradient of / must vanish at 0. 
Finally, observe that Re{/(#)} = — log \g(0)/z ( i\. By strict minimality of z, the modulus of g(w) = g{6) 
is greater than \zd\ for any w G T(z). □ 



We now prove Theorems ^lj, [3^ and |3.5| in reverse order. We see from Lemma [lj that proving any 
of these theorems amounts to evaluating the quantity 3 in equation ([|]^). From Lemma L2 we see that 
is a stationary point for the function / as long as r £ dir(z). The function / is in general complex 
valued, but we will see in Theorem 5.4 that it may be treated as if it were real valued, given the strict 
minimality of the zero guaranteed by Lemma 4.2 and the nonsingularity hypothesis. In particular the 
leading term of the integral in (^) is (27r)^ _1 ) //2 ^(0)r^ 1 d ^ 2 divided by the product of the square roots 
of the eigenvalues of the Hessian. Once we have identified ^(0) = *l>(0) as G(0)/(zrfi7rf), the theorem 
follows directly from Theorem 5.4 . 

Theorem |3.3| follows from the more explicit asymptotic development given in Corollary |5,3| . Finally, 



to prove Theorem 3.1, it remains to compute the quantity /"(0) in terms of the partial derivatives of H. 
First we compute the derivatives of g. 

Lemma 4.3. In a neighborhood of(z,w), if) and the derivatives of g are as follows. 



(4.6) 
(4.7) 



Hz 

1 



TT TT2 

TT _0—LH ' Z U 

n ZZ & jj n ZW I jjcy 12 WW 



G(z,w) 



-wH w (z,w) 



Proof. Differentiate the equation H(z,g(z)) = to get H z + g'(z)H w = which is the same as (| 
Differentiate again to get 

H zz + 2g'H zw + g"H w + (g') 2 H w 

and use ( [4.6D to eliminate g', giving ( |4.7j ). The formula for ip follows from the definitions of ip and of the 
partial derivative. □ 
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PROOF of Theorem [D] via direct computation: We know from Lemma fO[ that / vanishes to order 
at least two at 0. To compute /"(0), observe first that /" — log g is linear in 6, so /" = (log<?)". When 
Z = ze l8 , we have (d/dO) = iZ(d/dZ), so 



f" = iZ 



d_ 

dZ 



iZ 



dlogg 
dZ 



Z- 



d (Zg' 



dZ \ g 



Expanding this yields 
(4.8) 



^ g' + Zg" Z\g'f 

£ 1 K 



By our assumption, G does not vanish at (z,w), so as long as /"(0) 0, we may use Theorem O to 
conclude that the leading term asymptotic for a r , s is the k = 2, 1 = term of (|3.5[) . The term bo there is 
equal to 

G(z,w) I 2 
~ -wH w (z,w)]j f"(0)' 

G(z,w) I 2 



V>(0)r/(0) =^)V 2 //"(0) 



Thus from Theorem 3.3, 



-4(2,0) 



z- r w~ s 



2tt wH w {z,w)]j sf"(0)' 

Now evaluate this using the value A(2, 0) = yjix and equation (|4.8|) along with (4.6) and ([ 



to obtain 



1 _ r _ s G(z,w) (-wH w (z,w)) 2 
z w 



wH w (z,w) 



sQ 



where 



Q = (-wH w (z,w)ff"(0) = {-wH w (z,w)fz 
With the help of Lemma Ell we see (using g(z) = w) that 



3 -g'( z )-zg"(z) , z\g'{z)) 



Q 



-wH,., 



1 



-wH v 



-wH,, 



H 7 



> £ 



9{z) 



~H ZW -\- 



Hi 
H 2 



z 2 H 2 



+ 



2 H 2 

- LJ - in 



w 



evaluated at (z, w), which simplifies to the expression in Theorem 3.1. We see also that the nonvanishing 



hypotheses on Q is enough to guarantee /"(0) ^ 0, which finishes the proof of Theorem 34. 



□ 



5. Some oscillating integrals 



The oscillating integrals we require are integrals over a neighborhood of zero in M. d of the complex- valued 
integrand: 

exp(-A/(x))V>(x)<ix 



N 

where /(0) = 0, V/(0) = and Re{/} > 0. They are not difficult to compute, but since the standard 
references assume / is either real or purely imaginary, we sketch the development of these results. We 
mostly follow the exposition of (Stein ( |1993[) , adapting it to complex- valued phase functions and simplifying 
it to take advantage of the decay of the magnitude of the integrand in this case. 

We begin with one-dimensional results. Let denote the class of smooth functions with compact 
support. The following proposition is a well known consequence of Watson's Lemma (see, for example, 
Wonj ( gggg , Ch. 2, Theorem 1)). 
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Proposition 5.1. Let ip G C^°(M) and denote bj = ip^(0)/j\. Then as A — > oo, there is an asymptotic 
development 



/ exp(-Ax fc )V>(x) dx ~ V A + (k, l)bi\- {l+1)/k 



where, as in (3.1). 

A+(k,l) := k~ l T ' 



□ 

We extend this to more general one-sided integrals by a complex change of variables. Given any 
analytic, complex- valued function / on an interval [0, B], suppose that /(0) = 0, that /' / on (0,-B], 
and let k > 1 be the minimal so that f^ k \0) ^ 0. Let ip G Cq° vanishing to order Z > at 0. Denote 
Cj = /0')(0)/j! and = V (j) (0)/j!- The real part of Cfc is necessarily nonnegative. Define a function y on 
[0..B] by 

where the argument of c^ k is between — n/(2k) and n/(2k). The quantity /(ic) — c^x k is 0(x fc+1 ) near 
zero, so y is analytic near 0, and, in particular, is a diffeomorphism between [0, B] and a contour 7 from 
to some B*. Let F invert y. The derivatives of F at are easy to compute formally and the first j + 1 
starting from the depend only on the first j coefficients of / starting at c^. Define 

i/j* = (tpo F) • F' ; 

(5.1) b* = r U) (0)/j\. 

Theorem 5.2. Let f be analytic (complex-valued) on an interval [0, B]. Assume that /(0) = 0, that 
/' / on (0,-B], and Re{/} has a strict minimum at 0. Let k > 2 be minimal such that f^ k \0) 7^ 
and m be minimal so that the real part of f( m \0) does not vanish. Let ip G Cq° , let I be minimal such 
that V (0 (0) ± 0, and denote Cj := f {j) (0)/jl, bj := ^(0)/jl. Define b* as in {[5lQ. Then there is an 
asymptotic development 

,-B 00 

(5.2) / exp(-A/(x))^(x)dx~ V^ + (A:,j)&*A-^ +1 )/ fc . 

The constant in the 0(\~( N+1 ^ k ) term depends continuously (only) on the derivatives of f and ip up to 
(N + l)m/k - 1. 

Proof. Changing variables to y = f(x) 1 ^ k , the integral becomes 

exp(-\y k )ip*(y)dy; 



the curve 7 is the image of [0, B] under y, so 7 ; (0) = c\j k and 7 remains in the right half plane, strictly 
except at 0. For < N < M write ip* as Pm + y M+l RM, where Pm is a polynomial of degree M and 
Rm is bounded; this can be done since ip* may be approximated by a degree M polynomial to within 
0(y M+1 ) at 0. 
First, evaluate 

exp(-Xy k )P M (y)dy 



ASYMPTOTICS OF MULTIVARIATE SEQUENCES I 



17 



by moving the contour. Replace 7 by two line segments, the first of which goes along the positive real 
axis to some distance e and the second of which is strictly in the right half plane (we assumed Re{/} > 
except at 0). The integral along the second segment is exponentially small since the integrand is. Hence 
the combined contribution is the series ( |5.2| ) out to the j = M term. 
Next, bound 



/ exp(-Ay fe )y M+1 i? M (y) 



dy 



With C representing different constants in different lines, we now observe that on 7 we have Re{— y k } < 
— C\y\ m . Thus, parametrizing 7 by arc-length, an upper bound is given by 



poo 

/ exp(-ACi m )i 
Jo 



m\,M+l 



C\R M in{t))\dt. 



This is easily seen to be bounded above by CX~ ( M + 2 )/ m where C depends on the first M derivatives of 
/ and if). Choosing M > m(N + l)k — 1 we have a remainder term that is 0(\~( N+1 ^ k ), proving the 
theorem. □ 



The value of a two-sided integral follows as a corollary. 



Corollary 5.3. Assume the hypotheses of Theorem 5.1, with f now defined on an interval 
there is an asymptotic development 

r-B 

(5.3) 



-B,B\. Then 



exp(-A/(x))^(x)dx A(k,j)b*\- {j+1)/k . 
B U 



with A(k,j) given by ( pO| ) and (3.S). The bounds on the remainder terms each depend continuously on 
finitely many derivatives of f and 1^ on [—B,B]. 

Proof. The two-sided integral is the sum of two one-sided integrals on intervals [0, B] and [—B,0]. The 
integral over [—5,0] may be written as an integral over [0, B] of the function exp(— A/(— x))ip(— x) dx. 
With b^ still denoting the coefficients resulting form the application of Theorem [|]^ to the first integral, 
let b*[ denote the coefficients when Theorem 5^ is applied to the second integral. In order to add the two 
integrals, we write bj in terms of b*[ by means of the following routine computation. 
Let Cfc := c k e ia with c > and |a| < ir/2 and define the analytic quantity R so that 



y{x)= c k x k {l + R{x)) 1 



i/fc 



ce ia/k x(l + R{x)). 



If k is odd, then then the hypothesis Re{/} > implies that is purely imaginary. We have 



y[x) 



-c k x k {l + R{-x)) 1 



l/k 



ce- ia/k x(l + R(-x)) = -y{-x)e 



-2ia/k 



Writing n for the inverse function to y and 57 for the inverse function to y we then have 

i){x) = — r/(- 
Hence, letting C/[-] denote the coefficient of y l , 



-e 2ia/k x) 



bj 



Ci [ip(-fj(x)) ■ ff{x)] 



Ci 



2ia/k 



x)) ■ e 



2ia/k _ 



-e 2ia/k x) 



and thus 



(_iy e 2ia(l+l)/k b * _ 
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When k is even, the computation is similar but easier, resulting in 

b* = . 

Now observe that if k is odd, hence is purely imaginary, then e 2ia ( l + 1 )/ k = e ±*7r(i+i)/ft acC ord- 
ing to the sign of the argument of c&. Setting A(k,l) = (1 + (— ly)A+(k, I) if k is even and (1 + 

we recover the definition in (|3.2|) and (^) and prove the Corol- 
lary. □ 

Theorem 5.4. Let f be a smooth complex-valued function on a neighborhood of in M. d such that 
Re{/} > with equality only at 0. Suppose further that V/(0) = 0, and that the Hessian (matrix 
of second partials) of f has eigenvalues with positive real parts. Let H denote the Hessian determinant at 
0. Then for tp € , there is an asymptotic expansion 

f exp(-A/(x))^(x)dx~^^A-^/ 2 

3>l 

where I is the degree of vanishing of t/j atO. If I = then Cq = ip(0)(2ir) d / 2 Ti.^ 1 / 2 . The choice of square 

root is determined by Tt^ 1 ^ 2 = YYj=i ^ 2 where [ij are the eigenvalues of the Hessian and the principal 
square root is taken in each case. 

Proof. Let Q = =i Q.i,j z i z j be the quadratic form determined by the Hessian at the origin. Denote 
the eigenvalues of Q by {/ij : 1 < j < d} and note that each fij has nonnegative real part. 

Step 1: change coordinates to make / exactly equal to the quadratic form Q. Indeed since /(x) = 
Q(x)/2 + 0(|x| 3 ), and the Hessian is nondegenerate, there is a locally smooth change of variables {xj(z) : 
1 < J ' <■ d} such that /(z) = Q(x(z))/2 and the Jacobian at the origin is 1. 

Step 2: normalize by 7i 1 ^ 2 . For any quadratic form Q there is a linear change of variables y(x) such 
that Q(x) = Yl^iVj- The change of variables matrix P satisfies PP T = M(Q), the symmetric matrix 
representing Q. Changing variables to y introduces an integrating factor of detP which is a square root 
of H since M(Q) is just the Hessian. Let J\f' be the region of integration over which y varies when z 
varies over an appropriately small neighborhood of 0. 

Step 3: Expand tp into monomials. The function ip has now become tp, where ^(0) = "HT 1 ! 2 ip(0) 
and the sign of the square root will be chosen later. We may expand ip into monomials, using the same 
argument as in the proof of Theorem |5.2j to show the remainder term can be made Odyl^) for any N. It 
remains to evaluate the integral over the region of integration, N' of 

r. Ct 

/ exp(-A yV 2 )VKy)<iy 

when ip is a monomial. 

Step 4: move the region of integration to the real d-space. Let M" be the projection of J\f' onto 
M. d by setting the imaginary part to zero. We claim that changing the region of integration from J\f' 
to N" alters the integral by an amount rapidly decreasing in A. To show this, let Q be the region 
{Re{x} + z£lm{x} : x G M',t £ [0,1]}. The boundary of SI (as a manifold) is composed of N\M" 
(with opposite signs) together with S := {Re{x} + itlm{x} : x G dftf',t £ [0,1]}. For any d-form u, 
J n doj = J dn oj. When to = exp(— A Y2j=i ^jVj ^dyi A • • • A dyd is a holomorphic (i-form, we see that doj 
vanishes (being the sum of d/dlj terms) so that 



ijj = / u + u. 

JV' JJV" Js 
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We know that Re{]T\ fijy'j} is bounded away from on dM', and its minimal value on S lies on cW 7 , 
hence the integral over S decays exponentially. 

Step 5: evaluate the integral. Factoring y r exp(— A Y2j=i Vj) one-dimensional integrals and 
plugging into Proposition [5T] yields an asymptotic expansion whose leading term (when I = 0) is equal to 
(27r) d / 2 V>(0) Ti.~ 1/2 . When /(z) is the function Y$=i z j> tnen tne positive square root is taken. The choice 
of square root must be continuous in the analytic topology on functions having nondegenerate Hessians 
and having eigenvalues with positive real parts, and the only such choice is the product of the principal 
square roots of the eigenvalues of the Hessian. □ 

6. Classification of cases 

For purposes of classification some natural questions are: 

(i) what are all possible local geometries of minimal points of V? 

(ii) which of these can be handled by variants of the methods in this paper? 

(hi) are these sufficient to yield a good approximation to a r no matter what the direction, r/|r|, and no 
matter which generating function in the class, say, of functions meromorphic in a neighborhood of 
their domain of convergence? 

To make the last question more concrete, consider the simplest possible example, namely binomial coef- 
ficients, where F = 1/(1 — z — w) and V is a complex line. There are no singular points here, but how do 
we know that as (z,w) varies over minimal points of V, the direction dir(z,-u;) will cover all of MP 1 ? 



This question will be answered by Theorem 3^, but first we need to add some detail to the geometric 
discussion begun in Section It will be evident that quite a few cases need to be considered, some 
of which require new tools and some of which require only minor modifications. Accordingly, the results 
will appear in several papers, currently under preparation. In other words, a discussion of taxonomy will 
necessarily refer to results not yet published, and we will indicate to the best of our knowledge which ones 
are expected to be routine. 

Given a point z E V, we extend the definition of dir(z) to mean the set of limits of dir(y) as y — ► z 
along smooth points. When z is minimal, this is just the set of normals to support hyperplanes of logP at 
the point (log |zi|, . . . , log |<Zd|), so this is consistent with the old definition. As we will see shortly, dir(z) 
may be a (d — l)-dimensional subset of PJP^ 1 when z is a critical point of V. 



When H has a repeated factor, the residue computation in equation (^4) must be replaced by one 
involving the derivative. The remainder of the computation proceeds without a hitch as before. Details 
are given in Pemantle Sz Wilson ( 2000bD . For the remainder of the taxonomy, we assume H to be square- 



free. Toral smooth points may be handled by methods exactly the same as strictly minimal points. The 
inner integrand in (|4.3| ) will in this case have its maximal modulus on a set of dimension larger than zero. 
A modification of the necessary oscillating integral computation that works in this case is also given in 



Pemantle fc Wilson| (pOOObP 



If z G V is not smooth, all the first partials vanish. The expansion of H(x) near z is then a sum of 
terms of degrees 2 and higher. We call z a homogeneous point of degree k if this expansion contains terms 
(xj — Zj) k for each j = 1, . . . , d, and contains no terms of total degree less than k. 

Lemma 6.1. If z is a locally minimal point of V with nonzero coordinates, and F is meromorphic in a 
neighborhood of z then z is homogeneous. 

Proof. Passing to F(zix±, . . . , .z^d) if necessary, we may assume z = 1. Setting Xj = 1 for all but one 
index j, we cannot obtain the zero function (by minimality), and so some term in the expansion around 
1 is a pure power of (xj — 1), and we denote the minimal degree such term by Cj(xj — \) k i . If z is not 
a homogeneous point, then there is some j for which some monomial has total degree lower than kj. 
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Assume without loss of generality that j = d. The function F{ X <j X • • • • % x j y) then has a minimal degree 
pure y — 1 term co(y — l) fc , /c := k d , and some term c'(x — l) a (y — l) b with a + b < k. In other words, the 
Newton Polygon of F(x, ... , x, y) around (1, 1) has a support line passing through (0, k) with slope — p/q 
in lowest terms, and p > q. It is well known that we may describe the solutions y{x) of the equation 

F(l + x,... ,l + x,l + y) =0 

as follows. Write 

H := (y - l) k (co + c x {y - l)~*(x - If + c 2 (y - l)~^[x - 1) 2 « + • • • + c s (y - l)~ sp (x - l) s<1 ) 

for the polynomial collecting all the terms on this support line. Then for each q th root of unity, u>, and 
each root A of Yl c s-j^ = 0, there is a solution y = \ l l p x q l p (uj + o(l)) as x — > 0. A proof may be found 
Bricskorn & Knorrei (1986). 



m 



Varying x over the set \ir — arg(x)| < vr/4, we see that the solutions y(x) must sometimes be in this set 
as well. For those x, the points (1 + x, . . . , 1 + x, 1 + y) will be in V D D(l) \ T(l), violating minimality of 
1. By contradiction, we have shown that no monomial in the expansion around 1 has lower total degree 
than any pure power term, hence 1 is minimal. □ 

Continuing the taxonomy, suppose that z is a homogeneous point of V of degree k>2. We say that z is 
a multiple point if V is locally the union of k analytic surfaces. Algebraically, this means that the leading 
(order k) terms in the expansion of H near z factors into linear pieces. If the homogeneous point z is not 
a multiple point, we say it is a cone point. When d = 2 there are no cone points, since any homogeneous 
polynomial in 2 variables factors completely over C. 

Our understanding of cone points is not yet complete, but an analysis involving cone points is underway 
in pohn &: Pemantlc| ( pOOCj ). For multiple points, most of the story is given in pemantle fc Wilson| ( p000a|) . 
In particular, the following theorem is proved there. 



Theorem 6.2 (Pemantle &; Wilson (2000a)). Let z be an isolated, minimal, multiple point of V with 
multiplicity k. Let S C MP _1 be the set of outward normals to support hyperplanes to \ogD at the point 
(log \ z\\, . . . , log \ zd\)- Then there is an integer p > and a polynomial function (ft : S — » M such that the 
asymptotic expansion 

(6.1) a r ~ z- r 0(r) ^ C 3 {r d ) k -v/ 2 -i/ 2 

j 

holds uniformly as r varies over compact subsets of the interior of S. □ 



The extension to toral multiple points is given in Pemantle & Wilson ( 2000b|) . If a multiple point is 



not isolated or toral, then the degree of multiplicity, k, must be less than the dimension, d. This cannot 
happen of course when d = 2, but does happen when d > 3. The method for handling this case, toral 



or otherwise, is given in Pemantle & Wilson (2000b). That paper will also contain some subcases of the 



isolated multiple point case, namely when the sheets of V intersect non-transversely. 

Having more or less completed the taxonomy, we now discuss when we can guarantee that our methods 
yield asymptotics in all directions. 

Theorem 6.3. Let F = G/H = ^ a r ^ s z r w s be the quotient of analytic functions G, H : C 2 — > C. Suppose 
that the coefficients a T)S are all nonnegative, and that F(z, 0) and F(0, w) are not entire. Then for every 
direction a E MP 1 there is a minimal z € V with a G dir(z). 

Proof. Let (x, y) be any point on the boundary of logP. For u < e x and v < e y the power series for F is 
convergent at (u, v). As u f e x and v f e y therefore, F(u,v) is finite and increasing. On the other hand, 
the power series for F is not absolutely convergent on T(e x , e y ), since we know F to have some singularity 



ASYMPTOTICS OF MULTIVARIATE SEQUENCES I 



21 



on this torus. Hence F(u,v) f oo as (it, v) f (e^e 1 '). Since .F is meromorphic, it must have a pole at 
(e x ,e v ), hence (e x ,e v ) £ V and is a minimal point of V. As (x, y) varies over the boundary of logD, we 
let 7 C V denote the curve traced out by this minimal point. 



Pick any a £ 



The convex set logP has horizontal and vertical support hyperplanes (by non- 



entirety of F(z, 0) and F(0, it;)), and therefore has a support hyperplane normal to a; let (x, y) be a point 
of intersection of this support plane with logX>. We have just seen that z(a) := (e x , e y ) is a minimal point 
of V. If z is a smooth point of V then a £ dir(z): either z is finitely minimal, in which case Theo rem |3.3| 
applies, or it is toral, in which case the toral version of this theorem from Pemantle & Wilson (2000b) 
applies. 

Assume now that z is not a smooth point. By Lemma |6.1| , z is a homogeneous point, and since d = 2, 
z is a multiple point. Theorem 6.2 then shows that a £ dir(z) in this case as well. This finishes the 
proof. □ 



7. Further details and open questions 



The theorems in this and subsequent papers give estimates that are uniform away from the boundary of 
the domain in which they are valid. In order for all of these to be patched together so as to give estimates 
valid now matter how r — ► oo, one must determine the bandwidth around the boundary for which the 
boundary estimates on either side hold. For instance, suppose (z, w) is a multiple point of degree 2 and 
that dir(z, w) is the set of slopes between 1/2 and 2. It appears that the asymptotic estimate in Pemantle] 
|fc Wilson| (^OOOal) holding near the line {s = 2r} can be written so it is valid out to s = 2r + c-^/r. If the 
estimate for the region s/r > 2 + e can be widened so it holds to s = 2r + c^fr and a description given 
that is valid in the regime (s — 2r)/y / r — ► c, then the estimates will patch together completely. 

Another natural question is the universality of the method when the coefficients have mixed signs. We 



conjecture that Theorem |6.3| still holds, in the sense that for every direction there is point z £ V for which 
integration near z yields correct asymptotics. What we know is that z may no longer be minimal. For 
example, if G = 1 and 

H = (1 - (2/3)w - (l/3)z)(l + (l/3)w - (2/3)z) 

then the point (3/2,3/4) is not minimal but yields asymptotics in the diagonal direction; one sees this 
by integrating along a deformed torus rather than along T(3/2,3/4). In fact we conjecture that such a 
deformation always exists, but the topology seems not transparent enough to yield an easy proof. 

The class of algebraic functions is in some ways almost as nice as the set of rational functions, and nicer 
than the meromorphic functions. For one thing, an algebraic function is determined by a finite amount 
of data, and may thus easily be input into a symbolic math package. Gao & Richmond ( |1992| ) give an 
analysis of algebraic and logarithmic singularities, but sometimes the relevant singularities for algebraic 



functions are poles. For example, in Larsen and Lyons' analysis (Larsen & Lyons 1999) of merge times 
for coalescing particles, they find an algebraic function of the form 



F(z,w) 



X(z,w) 



w 



1-vT 



with x analytic. The branch of the square root is chosen so that at the origin the denominator is 2, not 
0. There is a branchline at z = 1, but for all directions in MP 1 , there is a smooth pole on the curve 
w = 1 + vl — z yielding asymptotics in the desired direction. It is natural to ask when this will happen, 
and how one can tell effectively. Some questions of effectiveness are addressed in Pemantle Sz Wilson| 
( |2000a ) and Pemantle & Wilson ( 2000b| ), but there is probably substantial room for improvements on an 
algorithmic level. 
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